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ABSTRACT 

Strong lensing has developed into an important astrophysical tool for probing 
both cosmology and galaxies (their structures, formations, and evolutions). Now 
several hundreds of strong lens systems produced by massive galaxies have been 
discovered, which may form well-defined samples useful for statistical analyses. 
To collect a relatively complete lens redshift data from various large system- 
atic surveys of gravitationally lensed quasars and check the possibility to use it 
as a future complementarity to other cosmological probes. We use the distri- 
bution of gravitationally-lensed image separations observed in the Cosmic Lens 
All-Sky Survey (CLASS), the PMN-NVSS Extragalactic Lens Survey (PANELS), 
the Sloan Digital Sky Survey (SDSS) and other surveys, considering a singular 
isothermal ellipsoid (SIE) model for galactic potentials as well as improved new 
measurements of the velocity dispersion function of galaxies based on the SDSS 
DR5 data and recent semi-analytical modeling of galaxy formation, to constrain 
two dark energy models (ACDM and constant w) under a flat universe assump- 
tion. We find that the current lens redshift data give a relatively weak constraint 
on the model parameters. However, by combing the redshift data with the bary- 
onic acoustic oscillation peak and the comic macrowave background data, we 
obtain more stringent results, which show that the flat A CDM model is still 
included at la. 

Subject headings: Gravitational lensing: strong - (Cosmology:) cosmological pa- 
rameters - (Cosmology:) dark energy 



1. Introduction 

The discovery of strong gravitational lensing in Q0957+561 (Walsh et al.1979) opened 
up the vast possibility to use strong lens systems in the study of cosmology and astrophysics. 
Up to now, strong lensing has developed into an important astrophysical tool for probing 
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. Now several hundreds of strong lens systems produced by massive galaxies have been 



discovered, but only ~ 90 galactic-scale strong lenses with known redshift of source and 
image separation can form well-defined samples useful for statistical analyses. 

These well-defined strong lenses are particularly useful not only for cons training the 



statistical prope rties of galaxies such a s optical region ve l ocity dispersions (I Chad 12005 



Chae et al.l 120061 ) and galaxy evolutions ( IChae &: Mad 120031 ; lOfek et al.ll2003l ). but also for 
constraining cosmological para meters such as the present matte r density Q m and the equa- 



tion of state of da rk energy w (|Chae 



statistical sample (IBrowne et al. 



2003 



2003 



Mitchell et all 120051 ) . For example, the CLASS 



Chad 120031 ) containing 13 lenses which s trictl y sat 



i sfy w ell-defined selection criteria was first extensively used by IChae et al.l (120021 ) and IChae 
(120031 ). who found fl m ~ 0.3 in a flat cosmology with non-evolving galaxy populations. 
Then this CLASS sample was reanalyzed with the velocity disp ersion function (V DF) of 
ear ly-type galaxie s deriv ed from the SDSS Data Release 1 (DR1 (jSheth et al.l 120031 )) galax- 



ies 



Mitchell et al.l (12005( 1 . 



In our paper, we summarize a complete lens redshift s ample from various imaging sur- 



veys including the Cosmic Lens ALL-Sky Survey ( CLAS S: (IBrowne et al.ll2003l : iMyers et al. 
20031)) .the JVAS, the Sloan Digital Sky Survey (SPSS: rtcjguri et ali hood)). the PMN-NVSS 
Extragalactic Lens Survey (PANELS; (jWinn et al.l 1200 ll )) and the Snapshot survey, which 
accumulates 29 galactic-scale l enses so far to for m a well-defined radio-selected lens sample. 
Newly measured J1620+1203 (IKayo et al.l 120091 ) from SDSS is also included. 



One emph asis we need to mak e is the measurement of the VDF of galaxies. IChad (120051 ) 
found that the ISheth et al.l (120031 ) VDF of early-type galaxies underestimated their abun- 
dance based on the Wilkinson Microwave Anisotropy Probe (WMAP) 1st year cosmolog y 
JSpergel et al. 2003 ). Then, with a new method of classifying galaxies ( Park Sz Choi 2005 ). 
Choi et al.l (bo07t ) made a new V DF measurement of early-type ga laxies based on the much 
larger SDSS Data Release 5 (DR5; lAdelman-McCarthv et al] (ho07h ). Recently, Ichael (bood ) 
have determined the VDF of the late-type population using the Tu l ly-Fis her relation and 
SIE galaxy model, which matches relatively well that of ISheth et al.l (120031 ) . More recently, 
Chad (120101 ) introduced a correction term for high velocity dispersions and used the Monte 
Carlo method to separately generate the early-type and late-type galaxi e s and to derive a 
total VDF for the entire population of galaxies (see also in Bernardi et al. ( 201ol )). However, 



-3 - 



the simulated data points for the total population of galaxies are not fitted well by the VDF 
of the morphologically-typed galaxy populations, whi ch might be due to the errors in the 
adopted correlations between luminosity and velocity ( jChad |2010| ). 



Moreover, strong lensing has also been extensively applied to constrain dark energy, one 
of the most important issues of the modern cosmology ever since the observations of type la 
supern ovae (SNe la) firs t indicating an accelerating expansion of the universe at the present 
stage (IRiess et al.l Il998l ). Among diverse dark energy models, the most simple candidate 
for the uniformly distributed dark energy is considered to be in the form of vacuum energy 
density or cosmological constant (A). However, the cosmological constant is always entangled 
with (i) fine tuning problem (present amount of the dark energy is so small compared with 
fundamental scale) and (ii) coincidence problem (dark energy density is comparable to critical 
density today). Alternatively there exist other choices, for example, an X- matter c omponent 
which is characterized by an equation of state p = wp, where — 1 < w < (jZhul 11998 



Peebles fc Ratral 120031 ) . The goal of this work is to use the lens redshift test combined 
with the revised VDF of all-type galaxies based on the SDSS DR5 data to constrain two 
cosmological models. 

This paper is organized as follows. In Section [2] we briefly describe the analysis method 
including assumptions about the lens population. We then present the lens redshift data from 
various surveys in Section [3j We further introduce two prevalent cosmologies and show the 
results of constraining cosmological parameters using MCMC method in Section HJ Finally, 
we conclude and make a discussion in Section EJ 



Lens redshift test 



The statistical lensing model used in this paper incorporates the (differential) lensing 
probabilities of specific image multiplicities for the multiply-imaged sources using the SIE 
lens model. The primary assumption we make is that the distribution of early-type galaxies 
in luminosity is given by the Schechter (1976) form 

L\ & { L_\dL 
L* J 

Considering a power-law relation between luminosity (L) and velocity dispersion (a), i.e. 
j- = , we can describe the distribution of galaxies in velocity dispersion in the form 

of the modified Schechter function 



<p(L)dL = 4>* 



— ] exp 



(1) 



dn 



{a) da 



a) flSheth et al.l 120031 : Mitchell et al.ll2005l ) 
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where 0* is the integrated number density of galaxies, <r* is the characteristic velocity dis- 
persion, with the following relations: a = (a + 1)7, — 7, and 0* = (fi*T(a + 1). 

Following Eq. ([2]), the particular differential probability that a source at redshift z s be 
multiply imaged by a distribution of galaxies at redshift z\ with a image separation A# can 
be written as: 



5p = 



d 2 p 1 dp 



dzd(A6) ' dz 
1 P 



2 r[(a + 4)/0\ A0» \A6 



A9 



a/2+1 



cxp 



Ad 
~A0l 



13/2 



(3) 



where A6* is the characteristic image separation given by 



A0* = A8vr 



D{z,z s 
D{0,z s 



a* 



(4) 



and D(z 1 ,z 2 ) is the angular-diameter distance between redshifts z\ and z 2 - On the hypoth- 
esis that galaxies are not biased toward oblate or prolate shape, we choose th e dynamical 
normalization factor A ~ 1 for the singular isothermal ellipsoid (SIE) model ( jChad 12003k 



20051 ). Though isothermal mass m odel would, in ge neral, be too simplified to do accu- 
rate modelling of individual lenses ( IChae et al.ll2002l ). it is accurate enough as first-order 



appr oximations to the mean properties of g a 



e.g.(Kochanek 19931: Mao &: Kochanek 



1994 



axies for the analyses of statistical lensing 



19971 : iFassnacht & Cohen! Il998l : iRusin & Kochanekll2005h ). 



Rix et al.1 11994 : iKochanekl 11990 : iKing et al. 



Following I Chad (120031 ). the likelihood C of the observed image separations for A^l 
multiply-imaged sources reads 



In C = OJi In 5pi 



1=1 



(5) 



(9) 



Here Spi is the particular differential probability for the /-th multiply-imaged source. 
(g = e, s) is the weight factor given to the early-type or the late-type populations, which 



satisfies oj\ < +0J^ 



1. If the lensing galaxy type is unknown, we use u l 
(g = e, s) (jChadl2003l ). Accordingly, a "x 2 "is defined as: 



(9) 



Sp\ 9) /[5p\ e) +Sp. 



X 



2 - -21n£. 



(6) 



Notice that the \ 2 here is free of the dimensionless Hubble constant h, which makes 
it possible to perform as an individual cosmological probe besides SNela, CMB, BAO etc. 
The 'best-fit' model parameters are determined by minimizing the \ 2 [Eq.®] and confidence 
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limits on the model parameters are obtained by using the usual A% 2 (= % 2 — Xmin) static, 
where Xmin * s the global minimum \ 2 f°r the best-fit parameters. Our cal culating method is 



based on the publicly available package CosmoMC (ILewis &: Bridld 120021 ) . 



On the side of the measurement o f VDF, the first d irect measurement of the VDF of 
early-type galaxies was the SDSS DR1 (ISheth et al.l 120031 ): 



DRl 



[(4.1 ±0.3) x 10 -i */i d Mpc" 
88.8 ± 17.7km s _1 , 
6.5 ±1.0, 1.93 ±0.22], 



(7) 



Then IChoi et al.l (120071 ) obtained a new VDF based on the much larger SDSS DR5: 

(0*, a*,a,/3) DR5 = [8.0 x 10- 3 /i 3 Mpc- 3 , 

161 ± 5km s -1 , 

2.32 ±0.10, 2.67 ±0.07]. (8) 

Obviously, the revised DR5 VDF, which has been proved to provide an e fficient way to 
constr ain dark energy models combined with gravitational lensing statistics (jZhu &: Mauro 
2008al ). is quite different from the DRl VDF in the characteristic velocity dispersion at la. 
While early-type galaxies dominate strong lensing, late-type galaxies cannot be neglected. 
However, for the late-type galaxy population, th e dire c t mea surement of the VDF is com- 
plicated by the significant rotations of the disks. I Chad (120061 ) estimated all the parameters 
of equation ([2]) for the late-type population using the Tully-Fisher relation and SIE galaxy 
model: 



&*,ac, /3)late 



[1.13 x 10- 1 /i 3 Mpc^ 3 , 
133 km s" 1 , 0.3, 2.91], 



(9) 



which matches relatively well that of Sheth et al. ( 2003 ) who determined o\j; late ' ) using a 
Tully-Fisher relation. 

Currently, it is found that s imple evolutions do not significantly affect the lensing statis- 



tics if all galaxies are early- type (IMao fc Kochaneklll994l : lRix et alll994l : lMitchell et al.ll2005 



Capelo fc Natarajanl 120071 ). Many previous studies on lensing statistics without evolutions 
of the velocity functio n have also got appealing results that agree w ith the galaxy number 
count s (llm et al.ll2002l ) and the redshift distribution of lens galaxies (IChaell2003l : lOfek et al. 
20031 ). However, in this paper, we consider galaxy evolutions both for ea rly-type and late 



type galaxies wi th a recent semi-analytical model of galaxy formation (IKang et al 
Chae et alJliooS ): 



2005 



fo(z) = 0*, o (l + z) Un - a,{z) = <7*, (1 + z) u \ 



(10) 
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with the best-fit parameters of iKang et al.l (120051 ) 



;-0.229, -0.01) for early-type 



and 



'1.24, —0.186) for late-type galaxies. 



3. Lens redshift data 

Large systematic surveys of gravitationally lensed quasars provide a large statistical 
lens sample appropriate for studying cosmology. In this section, we summarize a sample 
both from CLASS, SDSS observations and recent large-scale observations of galaxies, which 
will be used as the input for the statistical lensing model described in section [2J Two main 
sources of the lens redshift data are the Cosmic Lens All-Sky Survey and the Sloan Digital 
Sky Survey Quasar Lens Search. 

A. CLASS 



As the largest completed galactic mass scale gravitational lens search project, the Cosmic 
Lens All-Sky Survey, along with its predecessor project: the Jodrell-Bank-VLA Astrometric 
Survey (JVAS) has confirm ed 22 multiply- imaged systems out of a total of 16521 radio 
sources (IBrowne et al.ll2003l ). Out of the entire CLASS sample including the JVAS sources, 
a subsample of 8958 sources containing 13 multiply-imaged syst ems constitutes a st atistical 
sample that satisfies well-defined observational selection criteria (IBrowne et al.ll2003l ). which 
can be used for the statistical analysis of gravitational lensing in this paper. 

In this work we need galactic-scale strong lens samples that satisf y well-defined observa - 
tional select ion criteria. With the well-defined selection criteria from IBrowne et al.l ( 120031 ); 
Chad ( 120031 ). two out of the four multiply-imaged sources in the JVAS, 0414+054 with a too 
steep spectral index and 1030+074 with the fainter-to-brighter image flux-density ratio of two 
images less than 0.1, are excluded from the final CLASS statistical sample. Meanwhile, we 
stress that the measured lens redshifts, source redshifts, image separations, image multiplici- 
ties and the lensing galaxy types (if determined), as shown in table HJ are all needed through 
the likelihood function defined in Section [2] to constrain cosmological parameters. Under 
this criterion, for the CLASS sources, B0850+054, B0445+123, B0631+519, B1938+666 are 
clearly excluded due to their unknown source redshifts, while for B0739+366, J2004.1349, 
redshifts both for source and lens are unavailable 0. Moreover, the measured image separa- 



1 We discard lens systems that do not have measured lens or source redshifts, which may possibly cause 
biases. For many multiply-imaged sources without measured source redshifts, a possible strategy is to take 
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tions of 1359+154, 1608+656 and 2114+022 are not used because the observed angular sizes 
are due to multiple galaxies within their critical radii. 

B. SDSS 



The Sloan Digital Sky Survey Quasar Lens Search (SQLS; lOguri et al.l (120061)') is a 



photo metric and spectroscopic survey covering nearly a quarter of the entire sky (jYork et al. 



20001 ). We try to find suit able lens samples from the optical quasar catalog of the Sloan 
Digital Sky Survey (SDSS; lYork et al.l ( 120001 )). The first complete len s sample from Dat a 
Release 3 selected from 22,683 low-redshift (0.6 < z < 2.2) is provided in llnada et al.l ( 120081 ). 
It consists of 11 lensed quasars satisfying the following well-defined selection criteria: l)The 
image separation is between 1" and 20" with quasars brighter than i = 19.1; 2) The flux ratio 
of faint to bright images is greater than 10 -0 5 for double lenses. In this paper, an additional 
cut is applied to select a appropriate subsample: the lensing galaxy should be fainter than 
the quas ar components, becau se a too bright lens galaxy will strongly affect the colors of the 
quasars faichards et aJbooj ). Four lensed quasars, Q0957+561, SDSS J1004+4112, SDSS 
J1332+0347 and SDSS J1524+4409 are removed with this cut. Like the CLASS, due to the 
absence of lens redshifts, two more quasars, SDSS J1001+5027 and SDSS J1021+4913 are 
excluded. However, we successfully add four lens redshift data including SDSS J1620+1203, 
one of the ei ght newly discover ed and confirmed two-image lensed quasars by SDSS Quasar 
Lens Search (IKayo et al.ll2009l ). 



In table [TJ we summarize 29 Strongly-lensed Sources (redshifts both for sources and 
lenses, as well as the largest image separations and galaxy types) from the CLASS (JVAS), 
the SDSS, the PANELS and the Snapshot, which formulate a well-defined combined statis- 
tical sample for model constraints. 



4. DE models and constraint results 

In what follows, we choose two popular dark energy models and examine whether they 
are consistent with the lens redshift data listed above. 



1. Cosmological constant model. 



z a = 2, which is the mean source redshift for the multiply-imaged sources with measured source redshifts. 
However, in order to ensure the accuracy of constraint, we choose to abandon such a choice in our paper. 
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Table 1: Summary of Strongly- lensed Sources. The properties of the str ongly lensed systems 
from the Snapshot optic al survey and th e CL ASS (Browne et al.l 120031 ) and PANELS radio 
surveys are revised from iKochanek! (119961 ) and IChad (120031 . 120051 ) . Refere nces a r e the follow- 



in g: 1 - the CASTLE S website (http : //cfa -ww w.harvard.edu/cast les /); 2 -IChad (1 2003 



3- 
7- 
11 



OscozetaL 



1997 ): 8 



Eigenbrod et al. 



Qguri et all (120041): 15 



21 



(120071): 18 -iQguri et al. 



Lubin et al 



J2006ah: 



12 



(20001); 9 



Surpi fc Blandford 



(2008 



19 



(12003 ): 22 



Inada et al.l (12003b): 13 



Walsh et al. (119791): 16 -lYoung et al.l /I198lh : 17 



Inada et all ( 2003al): 10 



Kochanekl dl996h : 4 - E^g et al.l (119971): 5 - llnada_et al.l (I2OO5I) 6 - lEigenbrod et al 



Ofek et al. 



Eigenbrod et al. (2006b D; 14 



Fassnacht fc Cohenl ( 



19981): 20 



Kavo et all (120091 ): 23 - IWinn et al.l (120011 ) 



2005); 



(2007); 



( 20031 ): 



Morokuma et al. 



McKean et al. 



(120071 ): 







Source 


Lens 


Maximum 


Number 


Lensing 


Source 


Survey 


redshift 


redshift 


image 


of 


galaxy (-ies) 








(zi) 


separation (") 

sr \ / 


images 


type 


B0414+054 


CLASS 


2.64 


0.96 


2.03 


4 


earlv-tvpe (1, 2) 


B0712+472 1 " 


CLASS 


1.34 


0.41 


1.27 


4 


earlv-tvpe (1, 2) 


B1030+074 


JVAS 


1.535 


0.599 


1.56 


2 


2Gs (E+?) (1, 2) 


B1422+231 1 " 


JVAS 


3.62 


0.34 


1.28 


4 


earlv-tvpe fl, 2) 


J1632-0033 


PANELS 


3.42 


1 


1.47 


2 


early-type (1, 2) 


J1838-3427 


PANELS 


2.78 


0.36 


1.0 


2 


earlv-tvpe (1, 2) 


B1933+503 f 


CLASS 


2.62 


0.755 


1.17 


4 


earlv-tvpe (1, 2) 


Q0142-100 


Snapshot 


2.72 


0.49 


2.23 


2 


early-type (1, 3) 

J J sr V ) / 


PG1115+080 


Snapshot 


1.72 


0.31 


2.43 


4 


early-type (1, 3) 


B1938+666 


CLASS 


1.8 


0.88 


0.91 


4+2+R 


early-type (1, 2, 4) 


J0246-0825 


SDSS 


1.685 


0.723 


1.04 


2 


early-type (5, 6) 

J J sr \ i / 


SBS0909+523 


SDSS 


1.377 


0.83 


1.11 


2 


early-type (7, 8) 


J0924+0219 


SDSS 


1.523 


0.393 


1.78 


4 


early-type (9, 10, 11) 


J1226-0006 


SDSS 


1.125 


0.517 


1.24 


2 


early-type (12, 13) 


J1335+0118 


SDSS 


1.571 


0.440 


1.57 


2 


early-type (14, 13) 


Q0957+561 


SDSS 


1.413 


0.36 


6.17 


2 


early-type (15, 16) 


J1332+0347 


SDSS 


1.438 


0.191 


1.14 


2 


early-type (17) 


J1524+4409 


SDSS 


1.210 


0.310 


1.67 


2 


early-type (18) 


0712+472 


CLASS 


1.34 


0.41 


1.27 


4 


early-type (19) 


B1359+154 


CLASS 


3.24 


1 


1.67 


6 


3Gs (E+?+?) (1, 2) 


B2045+265 


CLASS 


4.3 


0.87 


1.91 


4 


2Gs (E+?) (2, 20) 


B1608+656 


CLASS 


1.39 


0.64 


2.08 


4 


2Gs (E+L) (2, 21) 


B0128+437 


CLASS 


3.12 


1.15 


0.55 


4 


unknown (6, 14) 


61152+199+ 


CLASS 


1.019 


0.439 


1.56 


2 


2Gs [?(E)a+?] (1, 2) 


Q1208+1011 


Snapshot 


3.80 


1.13 


0.48 


2 


unknown (5, 7, 8) 


J1620+1203 


SDSS 


1.158 


0.398 


2.765 


2 


unknown (22) 


B0218+357 


CLASS 


0.96 


0.68 


0.33 


2 


late- type (1, 2) 


B1600+434 


CLASS 


1.59 


0.41 


1.38 


2 


late- type (1, 2) 


J0134.0931 


PANELS 


2.23 


0.76 


0.68 


5+2 


2Gs (L?+L?) (1, 2, 23) 
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2. Dark energy with constant equation of state. 

Both of the two classes of models are currently viable candidates to explain the observed 
cosmic acceleration. Unless stated otherwise, throughout the paper we calculate the best fit 
values found, and vary the parameters within their 2a uncertainties for either class of model. 
Next, we shall outline the basic equations describing the evolution of the cosmic expansion 
in both dark energy models and calculate the best-fit parameters. 



4.1. Constraint on the standard cosmological model (ACDM) 

4-1.1. Constraint from the lens redshift data 

In the simplest scenario, the dark energy is simply a cosmological constant, A, i.e. a 
component with constant equation of state w = p/p = — 1. If flatness of the FRW metric is 
assumed, the Hubble parameter according to the Friedmann equation is: 



H 



a 



+ a« 



where Q m and Q\ parameterize the density of matter and cosmological constant, respectively. 
Moreover, in the zero- curvature case (Q = Q m +Qj^ = 1), this model has only one independent 
parameter: 9 = ^a- 



We plot the likelihood distribution function for this model in Fig. [TJ The best-fit value 
of the parameter is: ^a = 0.85^Q^g. It is obvious that the lens redshift data only give a 
relatively weak constraint on the model parameter f^A, though the universally recognized 
value of Q\ = 0.75 is still included at 68.3% CL (la). To make a comparison, it is necessary 
to refer to the previous results: the cu rrent best fit valu e from cosmological observations 
is: ^a — 0.73 ± 0.04 in the flat case (IDavis et al.l 120071). w hich is in relatively stringent 
accordance with our result. Moreover, iKomatsu et al.l (120091 ) gave the best-fit parameter: 
tt m = 0.274 for the flat A CDM model from the WMAP 5-year results with the BAO and 
SN Union data. We find that the constraint result from the lens redshift data is marginally 
consistent with the previous works above. 
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4-1.2. Joint analysis with BAO and CMB 



For the BAO data, the parameter A is used, which, for a flat universe can be expressed 



as 



A 



(H(z BAO )/H o y/z 



zbao 



ZBAO 



dz 



H(z)/H l 



2/3 



(12) 



where z BAO = 0.35 and AisA = 0.469 ± 0.017 from the SDSS feisenstein et all 120051 ) . 



For the CMB data, the shift parameter 1Z is used, which may provide an effective way 
to constrain the parameters of dark energy models due to its large redshift distribution. 
Derived from the CMB data, 1Z takes the form as 



K 



2 CMB 



dz 



H(z)/H > 



(13) 



where zcmb = 1090 (IKomatsu et al.l 120091) is the redshift of recombination and the 5-year 
WMAP data give K = 1.710 ± 0.019 fcomatsu et~al"1l2009h . 



In Fig. [2] and |3J we show the result by combining the lens redshift data with CMB and 



CMB+BAO, respectively. The best-fit parameter i s: Q\ 



0.78±g;g| and fi A 



n 7 pr+0.02 



relatively satisfactory result consistent with that of iLi et al.l ( 120101 ) . Compared with Fig. [TJ 
it is apparent that fl\ is more tightly constrained with the joint data sets. 



4.2. Dark energy with constant equation of state (w) 

4-2.1. Constraint from the lens redshift data 

If allowing for a deviation from the simple w = — 1, a component with a constant value 
for the equation of state could be introduced. The accelerated expansion can be obtained 
when w < —1/3. In a zero-curvature universe, the Hubble parameter for this generic dark 
energy component with density Vt x then becomes: 

Obviously, when flatness is assumed (Q = Q m + = 1), it is a two-parameter model with 
the model parameters: 9 = {Q x , w}. 

For the redshift data only, the best-fit values of the parameters are: Q x = 0.78, w = 
—1.15. Fig. H] shows the confidence limits in the Q x — w plane. On the one h and, we have 



w < —0.52 at 68.3% CL, which is quite different from the result of lChad (120061 ): on the other 
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Fig. 1. — The likelihood distribution function for the A CDM model constrained 
by the lens redshift data. 




Fig. 2. — The likelihood distribution function for the A CDM model constrained by the lens 
redshift data combined with CMB. 
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Fig. 3. — The likelihood distribution function for the A CDM model constrained by the lens 
redshift data combined with CMB and BAO. 




Fig. 4. — The w model constrained by the lens redshift data: likelihood contours at 68.3% 
and 95.4% CL in the fl x — w plane. 
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hand, the Einstein's cosmological constant (w — —1) is still favored within la error region. 
Therefore, it seems that the present lens reds hift data do not necessarily favor the phantom 
DE model with w < — 1 (jCaldwell 2002.2003 ). However, it is still interesting to see whether 
this remains so with future lager and better lens redshift data. 



4-2.2. Joint analysis with BAO and CMB 



In Fig. [5] and El we plot the likelihood contours with the joint data by combining the 
lens redshift data with CMB and CMB+BAO in the Q x — w plane. The best-fit parameters 



are: fl T 



0.80±°;^, w 



— 1.12±5:|a and fi s 



71+ - 07 w 

u - ' i -0.07> w 



-°- 78 -a34- Notice that both 



Vt x and w are more stringently constrained with the joint observational data. Meanwhile, 
the currently pref erred values of w i n this model still include the cosmological constant case: 
w = —1.01 ±0.15 (IDavis et al.ll2007l ). Therefore, when the equation of state does not depend 
on the redshift, the dark energy is consistent with a flat cosmological constant model within 
la error region. 



5. Conclusion and discussion 



Recent new obser vations, such as SNela, Wilkinson Microwave Anisotropy Probe (W MAP) 
( iKomatsu et al.l 120091 ) and Ba ryon Acoustic Oscil lations (BAO) (jPercival et al.l 120091 ). the 
time drift of subtended angles (IZhang fc Zhull2009l ). the updated Gamma-ray bursts (GRB) 



(e.g., iGao et al.l 12010c iLiang et al.l 120101 ) have provided many robust tools to study the dy- 
namical behavior of the universe. However, it is still important to use other different probes 
to set bounds on cosmological parameters. In this work, we have followed this direction and 
used the distribution of gravitationally-lensed image separations observed in the Cosmic Lens 
All-Sky Survey (CLASS), the PMN-NVSS Extragalactic Lens Survey (PANELS), the Sloan 
Digital Sky Survey (SDSS) and other surveys to constrain cosmological models with the new 
measurements of the velocity dispersion function of galaxies based on the SDSS DR5 data 
and recent semi-analytical modeling of galaxy formation. Two dark energy models (ACDM 
and constant w) are considered under a flat universe assumption. 

For the zero-curvature A CDM model, although the lens redshift data can not tightly 
constrain the model parameter Q\ = 0.85^11, a stringent constraint can be obtained by 
combining the lens redshift data with the comic macrowave background data: £l\ = 0.781q q3 
and the baryonic acoustic oscillation peak data as well: £7a = 0.75to'o2- Furthermore, we 
consider a flat cosmology with a constant w dark energy. For the lens redshift data, we 
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Fig. 5. — The w model constrained by the lens redshift data combined with CMB: likelihood 
contours at 68.3% and 95.4% CL in the Q x — w plane. 

-0.4 
-0.6 
-0.8 
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Fig. 6. — The w model constrained by the lens redshift data combined with CMB and BAO: 
likelihood contours at 68.3% and 95.4% CL in the fl x — w plane. 
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have w < -0.52 at 68.3% CL, a result different from that of IChad fl2006r ) with w < -1.2 
(68.3% CL), therefore, this strong lensing data do not necessarily favor a super- negative 
equation of state for dark energy. However, the Einsteins cosmo logical constant (w — — 1) is 
still included within la error region. Likewise, adding CMB and CMB+BAO does lead to 
further improvements in parameter constraints with Q x = 0.80 



+0.17 
-0.17; 



W 



-1.12±?;|| and Q a 



0.71 



+0.07 
-0.07' 



w 



n 70+0.22 

"U. ^8_o.34; 



respectively. Therefore, it indicates that the cosmological constant 
model is still the best o ne to explain this lens redshift data, a conclusion in accordance 
with the previous works (IDavis et al.l 120071 ) and th e results from the WMAP and the large- 
scale structures in the SDSS luminous red galaxies (jSpergel et al.ll2003l ; iTegmark et al.ll2004j ; 
Eisenstein et al.l 120051 ). 



However, we also notice that, firstly, the implementation of singular isothermal ellipsoid 
model (SIE) may be a source of systematic errors. For example, a lens ellipticity of 0.4 can 
lead to a difference of AQ m rs —0.05 co mpared with the sph erical case due to the variation 
of magnification bias and cross section (IHuterer et al.l 120051 ) . Secondly, as for the source of 
the lens redshift data, we simply discard lens systems that do not have measured lens or 
source redshifts in this paper. This could also possibly bring biases and will be considered in 
our future work. Thirdly, though the lens redshift test applied in this paper is free from the 
magnification bias arising from the uncertain source counts, it may also lose the statistical 
power of absolute lensing probabilities (or "lensing rates"). Lastly, though we have used the 
VDF of galaxies based on a much larger SDSS Data Release in the strong lensing statistics, 
the accuracy of measurements on relevant parameters of VDF may also make a difference. 
Hopefully, Large new samples of strong lenses will be expected to be obtain ed in future 



wide-fi eld imaging surveys such as the Large Synoptic Survey Telescope (LSST; llvezic et al. 



( 120001 )) etc and withi n a few decades ne xt generation observation tools such as the Square 
Kilometre Array (e.g.. iBlake et al.l ( 120041 )) will also improve the precision of lensing statistics 
by orders of magnitude. Therefore, the lens redshift test can play an important role in 
uncovering the physical processes of galaxy formation and universe evolution with much 
larger and better lens redshift data. 

Summarizing, from the above discussion we may safely arrive at a conclusion that the 
results from the observational lens redshift data are relatively agreeable and furthermore the 
lens redshift test can be seen as a future complementarity to other cosmological probes. 
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